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Ol Line Emission in the Quasar PG 1116+215 
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Abstract 

By observing the near-infrared spectrum of the quasar PG 1116+215 at z = 0.176 and combining with 
the HST/FOS spectrum, we obtained the relative strengths of three permitted Ol lines (A1304, A8446, 
and A11287) in a quasar for the first time. The photon flux ratios of the Ol lines of the quasar were 
compared with those previously measured in a Seyfert 1 and six narrow-line Seyfert Is. No significant 
differences were found in the O I line flux ratios between the quasar and the other Seyferts, suggesting that 
the gas density in the O I and Fe II line-emitting regions in the quasar is of the same order as those in 
low-luminosity AGNs. It was also found that the line width of Ol A11287 is significantly narrower than 
that of Lya, which is consistent with Ol and Fell emission occurring in the partly ionized regions at the 
outermost portion of the broad-line region where velocities are small. 

Key words: galaxies: active — galaxies: individual (PG 1116+215) — galaxies: quasars: emission 
lines — line: formation 



1. Introduction 

The permitted Ol lines at 1304, 8446, and 11287 A are 
common features in AGN (active galactic nuclei) and have 
been used to study the physical properties of BLR (broad- 
line region) clouds (e.g., Grandi 1980; Kwan & Krolik 
1981; Rudy et al. 1989; Laor et al. 1997). As shown in 
the partial Grotrian diagram of Ol atom in figure 1, the 
transition from the ground state 2p 3 P to the excited state 
3d 3 D° has an excitation energy whose wavelength is A = 
1025.77 A, falling within the Doppler core of Ly/3(1025.72 
A) for gas at 10 4 K. The excited state 3d 3 D° decays 
either back to the ground state, or by the emission of a 
A11287 photon, to the intermediate excited state 3p 3 P. 
The latter decays by a A8446 photon to the lower excited 
state 3s 3 S°, which finally decays to the ground state by 
a A1304 photon. If Ly/3 pumping is the dominant process 
of O I line formation, the photon flux ratio of these three 
lines should be 1:1:1, and the photon flux ratio, especially 
between O I A1304 and A8446, can be used as a reddening 
indicator (Rudy et al. 1989; Laor et al. 1997; Rodrfguez- 
Ardila et al. 2002b). It can also be expected that if the O I 
line emission is dominated by Ly/3 pumping, the O I line 
intensities should be used to measure the microturbulence 
parameter in BLR clouds, because Ly/3 pumping would 
become more efficient in a larger microturbulent velocity 

* Visiting Astronomer, Kitt Peak National Observatory, National 
Optical Astronomy Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc. 
(AURA) under cooperative agreement with the National Science 
Foundation, USA. 



field. 

Several mechanisms that could modify the 1:1:1 pho- 
ton flux ratio have been suggested, such as the collisional 
excitation of O I A8446 and the Balmer continuum absorp- 
tion of Ol A1304 (Grandi 1983). Rodriguez- Ardila et al. 
(2002b) found that in six of their seven low-luminosity 
AGNs, Ly/3 pumping is not the only mechanism respon- 
sible for the Ol line emission, and the 1:1:1 photon flux 
ratio is altered. Their sample consisted of six NLSls (nar- 
row line Seyfert Is) and one Seyfert 1. Before Rodriguez- 
Ardila et al. (2002b), I Zw 1, classified as the NLS1, is 
the only object where the three O I lines have been stud- 
ied (Rudy et al. 1989, 2000; Laor et al. 1997). It is thus 
not clear whether the properties of the O I line emission 
found in low-luminosity AGNs can be applied to quasars. 
The O I line emission in quasars would differ from those of 
NLSls, because it is generally considered that the gas den- 
sity in BLR clouds is higher in NLSls than in quasars, and 
thus the collisional processes are more dominant in NLSls 
(Baldwin et al. 1988, 1996; Laor et al. 1997; Wilkes et al. 
1999; Kuraszkiewicz et al. 2000). Investigating the Ol 
line formation mechanisms would also provide knowledge 
about the physical properties of the Fe II emitting region, 
since the Fe II and O I lines are considered to be emitted 
in the same portion of BLR clouds (Rodriguez- Ardila et 
al. 2002a). Since Fell emission is the strongest coolant 
in BLR clouds, it would give an important clue to better 
understand the radiative mechanisms in BLR gas. 

In order to study the O I line formation mechanism in 
quasars, we have started a program of taking NIR(near- 
infrared) spectra of quasars. Combining the NIR spectra 
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with the HST/FOS ultraviolet spectra, we intend to ana- 
lyze the three Ol lines, i.e., Ol A1304, Ol A8446, and Ol 
A11287. We present the ratios of the three 1 lines of the 
quasar PG 1116+215 at z=0.176. This is the first quasar 
in which the three O I lines have been analyzed. 

2. Observation and Reduction 

2.1. Near- Infrared Spectroscopy 

The NIR spectrum of PG 1116+215, covering 0.9 
- 1.8/Ltm, including the redshifted Ol A8446 and 
Ol A11287 lines, was obtained using the longslit of 
FLAMINGOS (Florida Multi-object Imaging Near-IR 
Grism Observational Spectrometer; Elston 1998 ) on the 
KPNO (Kitt Peak National Observatory) 2.1m telescope 
on 2005 February 28. The detector is the Hawaii II 
2048 x 2048 HgCdTe science-grade array, divided into four 
quadrants with 8 amplifiers each. A 4-pixel slit with a 
scale of 0"606 pixel" 1 was placed on the object in the 
north-south direction under a seeing of 1". The first order 
of the JH grism was used, which gave a spectral resolution 
of 430 km s _1 . The object was shifted along the slit by 
20" between exposures. The total integration time was 
2400s, which consisted of eight 300s exposures. J-band 
photometry of the object was also performed for photo- 
metric calibration. 

Data reduction was performed using IRAF (Image 
Reduction and Analysis Facility). 1 OH airglow lines were 
used to calibrate the wavelength scale. The flux scaling 
was twofold: (1) the A-type star SAO 81808 was observed 
to calibrate the relative response within the J-band, and 
(2) the standard star AS 20-0, having J = 9.55 ± 0.01 
mag (Hunt et al. 1998), was used to determine the J- 
band magnitude. A 9500K blackbody was fitted to the 
spectrum of SAO 81808 to obtain the sensitivity curve 
within the J-band. Both standard stars were observed at 
an airmass of 1.1, which is similar to that of the object. 
The total uncertainty in flux calibration is 5%. The final 
spectrum was re-binned into a 10 A step, corresponding 
to a spectral resolution of 300 km s _1 . 

2.2. Ultraviolet Spectroscopy 

The HST (Hubble Space Telescope) ultraviolet spec- 
trum of PG 1116+215 was observed on 1993 February 19 
and 20 by using the FOS (Faint Object Spectrograph) 
with three gratings (G130H, G190H, and G270H), which 
cover 1200 - 3200 A in the observed frame. The integra- 
tion times were 8724s for G130H, 1878s for G190H, and 
751s for G270H. 

The HST spectrum used in this work was obtained from 
Evans & Koratkar (2004). They recalibrated the raw 
archival spectra using the latest algorithms and calibra- 
tion data. Spectral data contaminated by intermittent 
noisy diodes and cosmic-ray events were identified man- 
ually and eliminated. They combined multiple observa- 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation, USA. 



tions of the same source, if available, in such a way that 
the resultant spectrum would have the highest possible 
signal-to-noise ratio. We removed prominent geocoronal 
and galactic absorption lines from the spectrum of PG 
1116+215 by interpolating from both sides of the feature. 
The spectrum of PG 1116+215 was re-binned into a 1 A 
step, corresponding to a spectral resolution of 300 km s" 1 . 
The total uncertainty in the flux calibration was less than 
5%. 

2.3. Galactic and Internal Extinction 

In order to correct the galactic reddening, we adopted 
Eb-v = 0.095 ± 0.015 mag from the extinction map of the 
Milky Way based on the far-infrared emission observed by 
IRAS and COBE/DIRBE (Schlegel, Finkbeiner, & Davis 
1998). The resolution of their extinction map is 6'1. De- 
reddening of our spectrum was performed by using a galac- 
tic extinction curve presented by Pei (1992). 

We also sought the effect of the intrinsic extinction of 
PG 1116+215. Popovic (2003) analyzed the HST/FOS 
spectrum of PG 1116+215, and found a flux ratio of 
Ha/H/3 = 2.94 ± 0.74. When corrected for galactic red- 
dening, as described above, the ratio became 2.70. Dong 
et al. (2005) showed that the flux ratio of the broad com- 
ponents of Ha to H/3 is Ha/H/3 = 2.97 ± 0.36 for their 94 
blue AGN samples. These AGN samples are considered 
to be free of intrinsic extinction. Judging from Ha/H/3 
of PG 1116+215, which is similar to those of these blue 
AGNs, the intrinsic extinction in PG 1116+215 appears 
to be little. 

2.4- Variability 

It is widely known that Seyferts and quasars are highly 
variable. Since the observation dates of the UV and NIR 
spectra are separated by more than 12 yrs, a possible vari- 
ability effect should be carefully examined before compar- 
ing the Ol lines in the UV and NIR wavelengths. J = 
13.59 ± 0.03 mag is found in the 2MASS (Two Micron 
All Sky Survey) database, 2 while J = 12.96 ± 0.03 mag 
was obtained for our observation. It is thus clear that PG 
1116+215 varied its brightness by about 0.6 mag since 
the 2MASS observation on 1998 February 2. However, 
because the O I lines are formed in the outermost portion 
of the BLR, it is expected that their variability is suffi- 
ciently small to be ignored, even if the power-law contin- 
uum, which is supposed to be direct light from the central 
source, shows significant variability. In fact, a monitoring 
observation of the well-known variable Seyfert NGC 5548 
revealed that when the continuum varied its brightness 
by about three times, the Mgll line, which is emitted in 
the same region as the Ol lines, varied by less than 6% 
(Dietrich & Kollatschny 1995). These arguments suggest 
that the variability has little effect on the O I line fluxes 

2 This publication makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the 
National Science Foundation, USA. 
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Fig. 1. Partial Grotrian diagram of O I. The solid lines are 
used for permitted transitions, while the dashed lines are for 
semi-forbidden transitions. 

of our data. 

3. Results and Analysis 

Figure 2 shows the UV spectrum of PG 1116+215 on the 
top panel and the NIR spectrum on the bottom. The three 
1 lines are clearly detected. Note that the signal-to-noise 
ratio of the NIR spectrum at the wavelength range shorter 
than lfim is low because of the aberration at the edge 
of the detector. The Ol A11287 line flux, which is free 
from blending with other lines, was measured by integrat- 
ing all counts above the estimated continuum level in the 
wavelength range between 1.319/im and 1.338^m in the 
observed frame. The line width was measured by fitting a 
single Gaussian to the observed line. The FWHM of 1 
A11287 is 2000 ± 200 km s _1 , and significantly narrower 
than that of Lya (~ 5500 km s^ 1 ). Rodriguez- Ardila et 
al. (2002a) showed that the line widths and profiles of 
Fell, Ol, and Call in their NLS1 samples are very simi- 
lar, while they are significantly narrower than that of Pa/3. 
Their results are consistent with the Fell, Ol, and Call 
lines occurring in the partly ionized regions, 3 which are 
farther from the central ionizing source and have smaller 
velocities. Our narrow Ol A11287 line width relative to 
Lya implies that this geometrical configuration found in 
the NLSls also applies to quasars. 

Although O I A8446 is blended with the Ca II lines in 
the red wing, the blueward part of the line is almost free 
from its contamination. We thus fitted a single Gaussian 
to the line in the blueward part where the contamination 
by Call is not severe, and then measured the flux within 

3 Note that O I, Fen, and Can emitting regions should be heav- 
ily overlapped because of the similar ionization potentials, i.e., 
13.6eV for O I, 16.2eV for Fen, and 11.9eV for Can. 
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Fig. 2. UV spectrum taken with the HST/FOS shown on the 
top panel and the NIR spectrum on the bottom panel in the 
observed frame. The positions of the three Ol lines, Can 
lines, and Si n A1264 line are indicated by vertical lines. 

the fitted Gaussian. 

The broad feature at 1304 A results from the blending 
of an O I triplet at A1302.17, A1304.86, and A1306.0 and a 
Sill doublet at Sill A1304 and A1309. The total blend 
flux of the A1304 feature was measured by integrating 
all counts above the estimated continuum level between 
1524 A and 1548 A in the observed frame. The Ol and 
Sill lines are severely blended and de-blending with high 
spectral resolution is hopeless for quasars. However, the 
Sill doublet blended in the 1304 A feature is expected 
to be accompanied by considerable Sill line emission at 
1264 A (Constantin et al. 2002). Actually, Sill A1264 is 
clearly detected in UV spectrum of PG 1116+215 (figure 
2) and its flux is 44 ± 5% of the blend flux of the 1304 A 
feature. The photoionization models by Kwan & Krolik 
(1979) and Netzer (1980) predict that the Sill doublet at 
1304 A and 1309 A has 16 - 28% of the flux of Sill A1264 
(Dumont & Mathez 1981), implying that 7 - 12% of the 
blend flux comes from the Si II doublet and 88 - 93% from 
the Ol triplet. Laor et al. (1997) found that 50 - 56% 
of the blend flux is due to the O I triplet in I Zw 1 and 
Rodriguez-Ardila et al. (2002b) found that the average 
portion of O I flux contribution to the blend flux is 75% in 
their three NLSls. In a sample of high-redshift quasars, 
Constantin et al. (2002) concluded that the A1304 fea- 
ture is due only to Ol lines, because Sill lines, which are 
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Table 1. Measured flux of O I A1304, A8446, and A11287 for 
PG 1116+215. 



Line 


Flux (10 14 erg s 1 cm 2 ) 


Ol A1304 


27.4±2.1 


1 A8446 


8.1±1.7 


Oi A11287 


4.6±0.6 



expected to be accompanied with Sill A1304 and A1309, 
are not seen. From the discussions presented above, we 
assumed that 90% of the A1304 of PG 1116+215 is due 
to Ol lines. We emphasize that assuming the portion of 
Ol lines as being 75% (Rodriguez- Ardila et al. 2002b) 
or 100% (Constantin et al. 2002) does not change our 
conclusions. 

We list the measured flux of Ol A1304, A8446, and 
A11287 in Table 1. 

4. Mechanism of O I Line Formation 

The photon flux ratios for PG 1116+215 of Ol 
A11287/OI A8446 (hereafter ROI ir ) and Ol A1304/OI 
A8446 (hereafter ROI uv ) are given in table 2 along with 
those for a Seyfert 1 and six NLSls. The O I A1304 flux of 
I Zw 1 was obtained from Laor et al. (1997), while A8446 
and A11287 were obtained from Rudy et al. (2000). ROI ir 
and ROI uv of other 5 NLSls and a Seyfert 1 were obtained 
from Rodriguez- Ardila et al. (2002b). Note that while the 
quoted errors for PG 1116+215 and I Zw 1 are la signif- 
icant, those for other samples are 2a significant. 

The average value of ROI ir for the NLSls is 0.73 ± 
0.22, where the la value, 0.22, reflects only the scatter 
from object to object, not including the errors measured 
for each object. ROI ir of PG 1116+215 falls within this la 
scatter of the NLS1 average. Even if we neglect Ton S180, 
which has a significantly larger value of ROIh than the 
other five NLSls, the average value of ROIh for the NLSls 
is 0.66 ± 0.15 and that of PG 1116+215 falls within the la 
scatter. Thus, ROL r of PG 1116+215 is not significantly 
different from those of other NLS1 samples. In the same 
way, we can see no significant difference between the ROL r 
of the Seyfert 1 galaxy NGC 863 and those of the NLSls. 

ROL r of the AGN was first measured for I Zw 1 by Rudy 
et al. (1989). They found that ROIi r of this object is 
equal to unity, and suggested that 1 lines are formed by 
Ly/3 pumping. However, Rodriguez-Ardila et al. (2002b) 
found that ROIh of I Zw 1 falls below unity, using the 
spectrum with higher resolution and wider spectral cov- 
erage published by Rudy et al. (2000). They also exam- 
ined the spectra of six NLSls including I Zw 1, and one 
Seyfert 1, to find that in six of their seven samples, ROLr 
are significantly below unity. Our result for PG 1116+215 
indicated that this trend is also present in quasars. They 
discussed that the collisional excitation to the upper level 
of the O I A8446 transition enhances the strength of the O I 
A8446 relative to that of the O I A11287, so that ROI ir falls 
below unity. This conclusion was drawn from the presence 
of Ol A7774 in 1H 1934-063, whose relative strength to 



O I A8446 is consistent with the combination of Ly/3 pump- 
ing and collisional excitation mechanisms. 

Rodriguez-Ardila et al. (2002b) ruled out continuum 
fluorescence as the additional mechanism from the lack 
or weakness of other Ol lines, such as A7002, A7254, and 
A13165, which should be present if continuum fluorescence 
is at work (Grandi 1980). Measuring the strengths of O I 
A6048, A7774, and A7990, which are very weak relative to 
that of A8446, they also suggested that the contribution 
to Ol A8446 by recombination is no larger than a few 
percents. 

The average value of ROI uv for the NLSls is 0.33 ± 0.30, 
adopting the upper limit value for 1H 1934—063. ROI uv 
of PG 1116+215 and NGC 863 fall within this la scatter 
of the NLSls average, which means that the ROI uv values 
of these two objects are not significantly different from 
those of other NLS1 samples. These values are signifi- 
cantly below unity. Kwan & Krolik (1981) showed that 
the intrinsic ROI uv is 0.76, rather than 1, due to the colli- 
sional excitation of O I A8446 and the Balmer continuum 
absorption of O I A1304. Grandi (1983) found that O I 
A1304 could also be destroyed by collisional de-excitation 
of the O I A1304 transition and de-excitation of the upper 
term of O I A1304 via O i] A1641 and A2324 line emission. 
He calculated the fraction of populations of 3s 3 5° that 
result in the observable A1304, and found that up to half 
of the A1304 photon could be converted to Ol] A1641 be- 
fore it leaves the emission-line cloud. However, Laor et al. 
(1997) found that no strong line is detected at 1641 A for 
I Zw 1, and that such a line cannot add more than ~30% 
to the observed A1304 flux. Such an investigation was 
possible for I Zw 1 because the Hell A1640 line, which is 
usually the dominant feature around 1641 A, is blueshifted 
by about 10 A in the I Zw 1 spectrum. Thus, they ruled 
out this mechanism as the main contributor to destroy the 
O I A1304 photon. Consequently, the major mechanism of 
O I A1304 destruction is the Balmer continuum absorption 
or collisional de-excitation, or both. 

5. Discussion 

Laor et al. (1997) showed that the relative strength 
of Cm] A1909 of I Zw 1 is significantly lower than that 
observed in the typical AGN. They argued that this sup- 
pression of Cm] A1909 may imply that the typical BLR 
density in I Zw 1 is about an order of magnitude larger 
than in the typical AGN. This suppression of Cm] A1909 
is also seen in other narrow-line quasars (Baldwin et al. 
1988, 1996; Wilkes et al. 1999; Kuraszkiewicz et al. 2000), 
suggesting that the high density of the BLR gas is common 
in this type of AGN. However, it is not clear whether these 
arguments are applicable to the region where Fe II and O I 
lines are emitted. Kuraszkiewicz et al. (2000) showed 
that while line flux ratios of Lya, C IV A1549, Si IV A1440, 
Cm] A1909, and Sim] A1892 observed in the NLSls can 
be explained by 10-times lower ionization parameter and 
a few-times (<10) higher densities than the normal AGN, 
the Mgll line strength cannot be explained by these pa- 
rameter values. Since the Fe II and O I lines are considered 
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Table 2. Measured ROI uv and ROI ir for PG 1116+215, one Seyfertl, and six NLSls. 



Object 


Type 


Ms 


R01 U v'' + 


KUli r '' + 


References 


PG 1116+215 


quasar 


-24.7 


0.52 ± 0.12 


0.76 ± 0.19 




NGC 863 


Seyfertl 


-20.9 


0.26 ± 0.04 


0.55 ± 0.08 


1 


1H 1934-063 


NLS1 


-17.1 


< 0.19 


0.64 ± 0.05 


1 


Ark 564 


NLS1 


-20.3 


0.07 ± 0.01 


0.82 ± 0.03 


1 


Mrk 335 


NLS1 


-21.0 


0.20 ± 0.02 


0.64 ± 0.05 


1 


Mrk 1044 


NLS1 


-19.5 


0.43 ± 0.07 


0.42 ± 0.05 


1 


Ton S180 


NLS1 


-22.6 


0.91 ± 0.15 


1.08 ± 0.16 


1 


I Zw 1* 


NLS1 


-22.7 


0.19 ± 0.02 


0.76 ± 0.11 


2, 3 



* These values were obtained by modifying the entries in the catalog presented by Veron-Cetty & Veron (2003), 
assuming cosmological constants as i/o=70km/s, Om=0.3, and H,a=0.7. Originally assumed values are -f/o=50km/s 
and qo=0. 

f ROIuv and ROI ir are defined as the photon flux ratio of A1304/8446 and A11287/8446. 

* The quoted errors are la significant for PG 1116+215 and I Zw 1, while those of other samples are 2a significant. 

* We assumed that 53% of the A1304 feature is due to 1 lines, which is the average portion of optically thick and 
thin cases presented by Laor et al. (1997). 

u References.— (1) Rodriguez- Ardila et al. (2002b). (2) Laor et al. (1997). (3) Rudy et al. (2000). 



to be formed in the same region as the Mgll line (see, e.g., 
Kwan & Krolik 1981), the physical properties in the Fell 
emitting region could be quite different from the Cm] 
emitting region. 

Three Ol lines presented here (A1304, A8446, and 
A11287) provide a good indicator of the gas density in the 
Fell emitting region. This is because these Ol line emis- 
sions are seriously affected by several mechanisms that 
are sensitive to the gas density, as described above. Our 
observation of these three O I lines in a quasar enables, 
for the first time, to investigate whether there is a signif- 
icant difference in the gas density of Fe II emitting cloud 
between NLSls and quasars. The results are quite differ- 
ent from the situation indicated by the Cm] A1909 line. 
As described above, there are no significant differences in 
ROL r and ROI uv among the NLSls, the Seyfert 1, and 
the quasar. This indicates that the physical properties of 
O I emitting cloud affect the O I line formation in the sim- 
ilar way in these three types of AGN. In other words, the 
collisional processes are working to a similar extent in the 
NLSls, the Seyfertl, and the quasar. Thus, our Ol ob- 
servations did not find any significant differences among 
these types of AGN in the gas density in the outermost 
portion of the BLR where the Fe II and O I lines are emit- 
ted. This result would provide some clues for modeling the 
environment of the Fell emitting cloud in NLSls, Seyfert 
Is, and quasars. 

NLSls are known as strong Fell emitters. The effi- 
ciency of Fe II emission is sensitive to several physical pa- 
rameters other than the gas density, in which the most 
sensitive parameters are the ionization parameter, micro- 
turbulence, and input spectral energy distribution (Netzer 
& Wills 1983; Wills, Netzer, & Wills 1985; Verner et al. 
1999, 2003). Our results suggest that the Fell emission 
enhancement in the NLSls may not be caused by the high 
density of the BLR gas. If this is the case, some of the 
other parameters should be quite different between NLSls 



and quasars. 

6. Summary 

We performed NIR spectroscopy of the quasar PG 
1116+215. By combining the NIR spectrum with the UV 
spectrum taken with the HST/FOS, we obtained three 
Ol lines (A1304, A8446, and A11287). We found that the 
line width of Ol A11287 is narrower than that of Lya, 
which is consistent with O I and Fe II emission occurring in 
the partly ionized regions at the outermost portion of the 
BLR. We also found that the photon flux ratio of the three 
Ol lines significantly deviate from 1:1:1, the expected ra- 
tio in the case of pure Ly/3 pumping formation. This 
strongly suggests the contribution of mechanisms other 
than Ly/3 pumping to the Ol line formation/destruction, 
for which the best candidates are the density-sensitive pro- 
cesses. Furthermore, the obtained photon flux ratio for 
PG 1116+215 is not significantly different from those of 
the NLSls and the Seyfert 1. This indicates that the gas 
density in Fe II and O I emitting regions are not signif- 
icantly different among NLSls, Seyfert Is, and quasars, 
although a larger gas density in the NLSls is indicated 
by the Cm] A1909 strengths. We also suggest that the 
physical parameters other than the gas density should be 
quite different between NLSls and quasars to account for 
NLSls as strong Fell emitters. 

We are grateful to the referee, Jack Baldwin, for useful 
comments to improve this manuscript. We thank the staff 
of KPNO for technical support and assistance with the 
observation. We would also like to thank H. Fukushi for 
her help. The trip of YM and SO to Tucson was financially 
supported by Research Center for the Early Universe, The 
University of Tokyo. 
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